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C^l ; Abstract 

We review the superfluid to Mott-insulator transition of cold atoms in optical lattices. The experimental signatures 
QQ ' of the transition are discussed and the RPA theory of the Bose-Hubbard model briefly described. We point out that 

, the critical behavior at the transition, as well as the prediction by the RPA theory of a gapped mode (besides the 

Bogoliubov sound mode) in the superfluid phase, are diflicult to understand from the Bogoliubov theory. On the 
other hand, these findings appear to be intimately connected to the non-trivial infrared behavior of the superfluid 
phase as recently studied within the non-perturbative renormalization group. 
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1. Introduction 

The realization of ultracold atomic gases has 
opened a new area of research in atomic physics. 
While early experiments have focused on quantum 
phenomena associated to coherent matter waves and 
superfluidity [T| (Bose-Einstein condensation, in- 



terference between condensates, atom lasers, quan- 
tized vortices and vortex lattices, etc.), the empha- 
sis has recently shifted towards strongly correlated 
systems '2 following three major experimental de- 
velopments: the realization of quasi-lD/2D atomic 
gases using strongly anisotropic confinement traps, 
the possibility to tune the interaction strength by 



Preprint submitted to Elsevier 



14 November 2008 



Feshbach resonances in the atom-atom scattering 
amplitude f3^ , and the generation of strong periodic 
potentials (analog to the crystalline lattice in solids) 
by optical standing waves. Thus quantum phenom- 
ena typical of condensed-matter physics have been 
observed in cold atomic gases 4 . 

The Mott transition is one of the most remark- 
able quantum phenomena due to strong correla- 
tions. Whereas in solids it corresponds to a metal- 
insulator transition driven by the Coulomb repul- 
sion, in bosonic cold atoms in optical lattices it is 
a transition between a superfluid (SF) and a Mott 
insulator (MI). This quantum phase transition has 
been observed in 3D, ID and 2D ultracold atomic 
gases I5I6I7I8] . 

2. The superfluid to Mott-insulator 
transition in cold atoms 

Ultracold neutral bosonic atoms can be stored in 
magnetic traps and cooled down to very low tem- 
peratures {'^ 10 nK) where they condense into a su- 
perfluid below a critical temperature Tc ~ 100 ^K. 
The importance of the atom-atom interactions can 
be estimated from the ratio between interaction and 
kinetic energies, 7 ~ n^^^a, which can be expressed 
in terms of the particle density n and the s-wave 
scattering length a. For Rb atoms in a magnetic trap 
(n ~ 10^* cm"'^, a ^ 5 nm), 7 is typically of or- 
der 0.02, which corresponds to the weakly interact- 
ing dilute limit. The physical properties of the gas 
(collective modes, vortices, etc.) are well described 
by a macroscopic wave function satisfying the well- 
known Gross-Pitaevskii equation ^ . 

The importance of interactions can be drastically 
enhanced by subjecting the atomic gas to a periodic 
lattice potential, which can be created by counter- 
propagating laser beams. The interference between 
the two laser beams forms an optical standing wave 
whose electric field induces a dipole moment in the 
atom and leads to an interaction energy VbL(i") = 
— iQ;(a;L)|E(r)p with E(r) the electric field at posi- 
tion r (a(wL) denotes the polarizability of an atom). 
By using different arrangement of standing waves, 
one can create various optical lattices. Three orthog- 
onal standing waves correspond to a 3D cubic lattice 
and an optical potential 

VoUr) = Vo[siii^{kx) + sin2(fcy) + sin^ {kz)], (1) 

where k — 27r/A is the wavevcctor of the laser light. 
Vq represents the lattice potential depth and is di- 
rectly related to the laser light intensity. Provided 
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Fig. 1. Mean-field phase diagram of the Bose-Hubbard model 
showing the SF phase and the MI phases at commensurate 
filling n. The dashed lines correspond to a fixed density 
n = 0.2, n = 1 and n = 2. 2 denotes the number of nearest 
neighbors. 

the atoms remain in the lowest Bloch band of the 
lattice (which requires Vq to be large enough wrt 
the single atom recoil energy Er = h^k'^ /2m), the 
interacting boson gas can be described by the Bose- 
Hubbard model [lOlllj defined by the Hamiltonian 

-^if^-£i)nt + —^n,{'<T't-'^), (2) 

i i 

where i/jj-^'* is a creation/annihilation operator de- 
fined at the lattice site i, fii — ij^jipi, and /i denotes 
the chemical potential fixing the average density of 
bosons. The on-site energy originates form the 
magnetic trap that confines the atomic cloud as well 
as the Gaussian shape of the laser beams. The hop- 
ping amplitude t and the local repulsion U are given 
by [12i 

so that the ratio t/U that controls the physics of the 
Bose-Hubbard model can be tuned by varying Vq, 
i.e. the intensity of the laser beams. 

In the absence of the confining trap (e^ = const), 
the phase diagram of the Bose-Hubbard model can 
be obtained from simple arguments. When the ki- 
netic energy dominates {t/U ^ 1), the ground state 
is a superfluid. In the opposite limit of a strong 
lattice potential {t/U ^ 1), the interaction energy 
dominates and the ground state is a Mott insula- 
tor when the average number of atoms per site is 
integer. For non-integer fillings, the ground state is 
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always a superfluid irrespective of the strength of 
the interaction lOJ (Fig.[T]). The situation is a httle 
more subtle when the confining trap is taken into 
account [H]. In a local density approximation, the 
trap potential can be accounted for by a spatially 
varying chemical potential fj,i = ^ — Ci (taking = 
at the center of the trap). When moving from the 
trap center to the edge, the local chemical potential 
decreases from /i to zero. We then expect to observe 
all phases which exist for a chemical potential below 
fj, (and the same value of t/U) in the absence of the 
trap. If, for instance, the chemical potential /i falls 
into the n = 2 lob in Fig. [U the ground state of the 
gas will correspond to a n = 2 Mott region at the 
center of the trap surrounded by a SF region with 
spatially varying density l<n<2,an=l Mott 
region and, near the boundary of the trap, a SF re- 
gion with density n < 1. The fact that the density 
remains constant in the MI regions, while the local 
chemical potential varies, is a consequence of the in- 
compressibility of the Mott phase {dn/d(x = 0). 

The SF-MI transition in cold atoms was first ob- 
served by Greiner et al. in a 3D ^^Rb gas [5'6]. It 
has since then also been observed in ID and 2D 
gases 7 8. The transition is detected by switching 
off simultaneously the magnetic (confining) and op- 
tical potentials and allowing for typically t = 10 — 20 
ms of free expansion. The density distribution 7i(r, t) 
of the expanding cloud is then obtained by absorp- 
tion imaging. It can be expressed as [13l6j 



|i&(k)pTik 



k=.; 



(5) 



(6) 



where u;(k) is the Fourier transform of the Wannicr 
function w(r) corresponding to the lowest Bloch 
band of the optical lattice. n(r, t) therefore provides 
a direct measure of the momentum distribution 
'^k = (V'Uk)''/'(k)). In the superfluid state with a 
finite condensate fraction /, the momentum dis- 
tribution rik exhibits sharp peaks - analogous to 
the Bragg peaks in the static structure factor of a 
solid - at the reciprocal lattice wavevectors G. (In 
a homogeneous superfluid, / = no/n with no — 
lim|r|_oo(V'iV'o) = |(V'(k = 0))p/F the condensate 
density.) In the MI where {tpltpo) decays exponen- 
tially over a few lattice spacings, the "Bragg peaks" 
are significantly suppressed and broadened [12] . 
Thus, the superfluid phase is signaled by a high- 
contrast interference pattern, as expected for a 



periodic array of coherent matter- wave sources [S]. 
Recent measurements of the momentum distribu- 
tion function rik in the Mott phase of a 2D atomic 
gas [S] (see also [5]) are in excellent agreement (with 
no adjustable parameter) with the RPA theory [T3] 
discussed in the next section. From the momentum 
distribution profile, it is also possible to extract 
a so-called "coherent fraction" representing the 
weight of the sharp peaks and closely related to the 
condensate fraction /, and thus locate the position 
of the SF-MI transition [8J. 

Using an experimental technique based on spa- 
tially selective microwave transitions and spin- 
changing collisions. Foiling et al. have directly ob- 
served the formation of the spatial shell structure 
in the SF-MI transition [15116] . This technique en- 
ables to record the spatial distribution of lattice 
sites with different filling factors. As the system 
evolves from a superfluid to a Mott insulator, it 
reveals the formation of a distinct shell structure, 
in agreement with the qualitative discussion given 
above, and therefore definitively shows the existence 
of incompressible (Mott) phases. 

Another trademark of the Mott insulator is 
the existence of a gap in the excitation spectrum 
whereas the superfiuid phase is characterized by 
a gapless (Bogoliubov) sound mode. Deep in the 
MI the gap is given by U, and should vanish at 
the MI-SF transition. In the early experiment of 
Greiner et al. [5] , it was shown that the response of 
the MI to a potential gradient is compatible with 
the expected gapped spectrum. In principle, the ex- 
citation spectrum can be measured by two-photon 
Bragg spectroscopy. This technique has allowed to 
observe the gapless mode of a SF atomic gas [17], 
but for a gas in an optical lattice no convincing 
result has been obtained so far 1 71. 



3. RPA theory... and beyond 

Several theoretical approaches have been used 
to study the Bose-Hubbard model: mean-field 
theory |10I18I19] . random- phase-approximation 
(RPA) [14|20|21|22] . strong-coupling expansion [23], 
numerical calculations (Quantum Monte Carlo) 24] 
or variational approach [25] . In this section, we dis- 
cuss the RPA theory of the SF-MI transition as well 
as some open questions related to our understand- 
ing of superfluidity. 
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3.1. RPA theory 

The RPA is based on the mean-field decoupling 
of the hopping term in ([2]), 

^li'j ^ + r^i', - (7) 

(from now on we ignore the confining trap), where 
4'i — Taking (pi = (pQ as a uniform order pa- 

rameter, the Hamiltonian ([2|) becomes a sum of de- 
coupled one-site Hamiltonians which can be solved 
exactly |26j . A finite value of the order parameter 
(00 0) signals superfluidity. The resulting phase 
diagram in the {U /t, /i/i) plane in shown in Fig. [1] 
It is possible to go beyond this mean-field approx- 
imation in the following way. In the presence of a 
(fictitious) external source j* that couples to the bo- 
son operator ipi , the (imaginary-time) action reads 

•Srpa = •Sloe - / dT^^^2 ['^iii,j4>j + (l^lUjipj 

i 

with Sloe the local part (t = 0) and /3 = l/T [27l . To 
eliminate the dependence of the partition function 
Z[j*,j] on the external source, it is convenient to 
perform a Legendre transform to obtain the Gibbs 
free energy 

r[0*, 0] = rioc[r , 0] - / dTY,<l>:u.jc^j, (9) 

where rioc[0*,0] is the local part {t — 0). The 
value of the order parameter (po = (V'i)j*=j=o fol- 
lows from the equation of state (5r/5(/)i(r)|0Q = 
5r/6(f)*{T)\^g = and determines the condensate 
density hq = \(l)o\'^- The Green function is ob- 
tained from the second-order functional derivative 
of dH) [28 . In the RPA, it takes the simple form 

G-i(q,zc^„;</)o) = Gi;i(zc^„;0o)+t(q), (10) 

where t(q) is the Fourier transform of the hopping 
amplitude tij. Equation is typical of a strong- 
coupling expansion in t/U and becomes exact in 
the limit t -> [2"3TT4] . The local Green function 
Gioc(*'^n; 0o) should be calculated in the presence 
of an external source j* ,j such that — (po- Note 
that in the SF phase G and Go are 2x2 matrices 
with both normal (e.g. (tpip*)) and anomalous (e.g. 
{iptp)) components. 

Equation (fTO)) is the central result of the RPA the- 
ory of the Bose-Hubbard model. Superfluidity sets 
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Fig. 2. Excitation spectrum in 2D. Top: n = 1 MI {U/t = 30); 
bottom: SF phase (U/t = 20) [M]. The dotted hnes show the 
result obtained from the BogoUubov theory (which predicts 
the phase to be superfiuid). [F = (0,0), M = (7r,7r) and 
X = (7r,0).] 

in when the boson propagator becomes gapless, i.e. 
^ioc(*^" = 0) -I- t(q = 0) = 0. The resulting mean- 
field phase diagram is qualitatively correct in 3D and 
2D (Fig.[T|). Equation pO|) can also be used to obtain 
the momentum distribution rik fl4'. In the MI, the 
agreement with the experimental results is remark- 
able |6|8j . The excitation spectrum is obtained from 
the poles of the Green function G(q, a; -I- «O+;0o)- 
In the MI, one finds two gapped modes as shown in 
Fig. [2] for a 2D system. In the SF phase near the 
MI-SF transition, there are four excitation modes, 
two of which being gapless (sound modes) for q — > 
114120121122125] . In the limit f/ ^ and for a large 
number of bosons per site, the two gapped modes 
disappear and the RPA theory reproduces the re- 
sult of the Bogoliubov theory [22]. These findings 
suggest that the Bogoliubov theory might not be 
appropriate, even on a qualitative level, to describe 
the superfiuid phase near the SF-MI transition. 

3.2. Critical behavior at the SF-MI transition 

The limitations of the Bogoliubov theory are also 
apparent when one considers the critical behavior 
at the SF-MI transition. The critical theory can be 
studied from the action jl4|26j 
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J'^ dr J d'W[Z\V'<p\'^ + Zi^*drip 

+ y|c',VI'+r-o|V'l' + "l^l1- (11) 

The coefficient Zi is related to the inverse of the 
slope of the transition line in Fig. [T] and vanishes at 
the tip of the Mott lob. The MI-SF transition is then 
in the universality class of the XY model in d + 1 
dimensions with a dynamic critical exponent z — 1. 
It occurs at a fixed density and is accompanied by a 
vanishing of the excitation gap of the MI. Away from 
the tip of the Mott lob, Zi is nonzero and the tran- 
sition has a dynamic critical exponent z ~ 2 [TU]. 
The density is not conserved and only one of the two 
modes of the MI becomes gapless (i.e. the gap of the 
MI does not vanish at the transition). Not only does 
the Bogoliubov theory always predict the system to 
be superfluid (irrespective of the strength of the in- 
teractions), but it also corresponds to Zi = 1 and 
V = Q and therefore appears to be a rather poor 
starting point to understand the SF-MI behavior. 

3.3. Infrared behavior in the superfluid phase 

The Bogoliubov theory provides a microscopic ex- 
planation of superfluidity and many of its predic- 
tions have been confirmed in ultracold atomic gases. 
Nevertheless a clear understanding of the infrared 
behavior of interacting bosons at zero temperature 
has remained a challenging theoretical issue until 
very recently. Early attempts to go beyond the Bo- 
goliubov theory have revealed a singular perturba- 
tion expansion plagued by infrared divergences due 
to the presence of the Bose- Einstein condensate and 
the Goldstone mode [30]. These divergences cancel 
in most physical quantities but lead to a vanish- 
ing of the anomalous self-energy Y,i2{q) in the limit 
q = (q, w) — > although the linear spectrum and 
therefore the superfluidity are preserved [31]. This 
observation seriously calls into question the validity 
of the Bogoliubov theory where the linear spectrum 
relies on a finite value of E12 (£12(9) — A-Kano/ni). 
The physical origin of the vanishing of the anoma- 
lous self-energy is the divergence of the longitudi- 
nal correlation function which is driven by the gap- 
less (transverse) Goldstone mode ~ a general phe- 
nomenon in systems with a continuous broken sym- 
metry [32!. The coupling between longitudinal and 
transverse fluctuations is not taken into account in 
Gaussian fluctuation theories such as the Bogoli- 
ubov theory. 



The infrared behavior of zero-temperature Bose 
systems is now well understood in the framework 
of the non-perturbative renormalization group 
(NPRG) pIM] (see also Ref. ,35i). The NPRG ap- 
proach is based on an exact flow equation satisfled 
by the Gibbs free energy T[(f)*,(f)] (see Sec. 13. 1|) as 
fluctuations are gradually integrated out beyond 
the Bogoliubov theory. It reveals that the Bogoli- 
ubov form of the action is essentially modifled in the 
RG process when the spatial dimension d < 3. In 
the RG language, this means that the Bogoliubov 
flxed point is unstable when d < 3. The new flxed 
point is characterized by a "relativistic" action, i.e. 
an action of the type ([TT]) with Zi = and V 
whereas the Bogoliubov fixed point corresponds to 
Zi > and V = 0. (In practice the Bogoliubov the- 
ory remains valid in 3D systems as the RG flow is 
only logarithmic so that the new fixed point is not 
accessible in any finite size system.) 

It is quite remarkable that a proper treatment of a 
Bose gas in a continuum model (i.e. with no underly- 
ing periodic lattice) yields a low-energy action sim- 
ilar to that obtained from the Bose-Hubbard model 
near the MI-SF transition. In other words, the crit- 
ical behavior at the MI-SF transition appears to be 
intimately connected to the non-trivial infrared be- 
havior of the superfluid phase. A generalization of 
the NPRG technique to the lattice case should shed 
light on this issue. 

4. Conclusion 

Ultracold atomic gases allow us to study strongly- 
correlated systems in an unprecedentedly controlled 
manner. Not only do we have a detailed microscopic 
understanding of the Hamiltonian of the system 
realized in the laboratory, but the microscopic pa- 
rameters that control the physics can be tuned 
by varying external flclds. Hamiltonians typical of 
strongly-correlated systems (e.g. (Bose-)Hubbard 
models) can be simulated in cold atomic gases, and 
new systems - with no equivalent in condensed- 
matter physics - can also be created !T_ . 

A great success in the study of cold atoms has 
been the observation of the SF-MI transition. In this 
paper, we have reviewed the main experimental sig- 
natures of the SF and MI phases. The RPA theory 
of the Bose-Hubbard model provides a qualitative 
and sometimes quantitative description of the sys- 
tem. What has not been observed so far is the excita- 
tion spectrum and the critical behavior at the tran- 



5 



sition. Whether the experimental difhculties and the 
finite size (as well as the non-uniform density) of 
the atomic clouds will allow these observations is an 
open question. 

From a more theoretical side, we have pointed 
out the difficulty to understand the superfluid phase 
near the SF-MI transition, as well as the critical be- 
havior at the transition, from the Bogoliubov the- 
ory. On the other hand, the critical behavior ap- 
pears to be intimately connected to the non-trivial 
infrared behavior of the superfluid phase as recently 
studied within the non-perturbative renormaliza- 
tion group |33I34| . 
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